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Abstract. The semi-discrete convolution with the Box Sphne is an 
important tool in approximation theory. We give a formula for the 
difference between semi-discrete convolution and convolution with the 
Box Spline. This formula involves multiple Bernoulli polynomials . 

Key words: polynomial interpolation, box splines, zonotopes, hyper- 
plane arrangements, Bernoulli polynomials. 

1. Box Splines and semi-discrete convolution 

Let y be a n-dimensional real vector space equipped with a lattice A. If 
we choose a basis of the lattice A, then we may identify V with and A 
with Z". We choose here the Lebesgue measure dv associated to the lattice 
A. 

Let X = [ai,a2, ■ ■ ■ ,a]\f] be a sequence (a multiset) of non zero vectors 
in A. 

The zonotope Z{X) associated with X is the polytope 

N 

Z(X) := ^hai;ti € [0,1]}. 

i=\ 

In other words, Z{X') is the Minkowski sum of the segments [0, Cj] over all 
vectors ai € X. 

We denote by C[V] the space of (complex valued) polynomial functions 
on V. 

Recall that the Box Spline B{X) is the distribution on V such that, for 
a test function test on V, we have the equality 



(1) {B{X),test) = [ ■■■ I test(S2tiai)dti---dtN- 

Jti=0 JtM=0 j=l 

We also note {B{X),test) = fy B{X){v)test{v). 

The distribution B{X) is a probability measure supported on the zono- 
tope Z{X). If X is empty, then B{X) is the 5 distribution on V . For the 
basic properties of the Box Spline, we refer to [5] (or [6], chapter 16) . 

If D is any distribution on V , the convolution B{X) * D is well defined 
and is again a distribution on V. If D = f{v)dv is a smooth density, then 
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Figure 1. Affine topes for X = [ei, 62, ei + 62] 
BijC) * D = F{v)dv is a smooth density with 

F{v)= ••• / fiv-y2ii('i)dti---dtN. 

If X generates V, the zonotope is a full dimensional polytope, and B(X) 
is given by integration against a locally -function. Let us describe more 
precisely where this function is smooth. 

We continue to assume that X generates V. An hyperplane of V generated 
by a subsequence of elements of X is called admissible. An element of V is 
called (affine) regular, if no translate v + \ oi v hy any A in the lattice A 
lies in an admissible hyperplane. We denote by V5-eg,aff the open subset of 
V consisting of affine regular elements: the set T^eg.aff is the complement of 
the union of all the translates by A of admissible hyperplanes. A connected 
component r of the set of regular elements will be called a (affine) tope (see 
Figure [T|) . 

The choice of the Lebesgue measure dv on V allows us to identify distribu- 
tions and generalized functions: if F is a generalized function, Fdv is a dis- 
tribution. If the distribution Fdv is given by {Fdv, test) = Jy f{v)test{v)dv, 
with f{v) locally L^, we say that F is locally L^, and we use the same no- 
tation for F and the locally function /. 

A generalized function 6 on y will be called piecewise polynomial (relative 
to X, A) if: 

• the function b is locally L^, 

• on each tope r, there exists a polynomial function 6(r) on V such that 
the restriction of 6 to r coincides with the restriction of the polynomial 6(r) 
to r. 

If F is a piecewise polynomial function, we will say that the distribution 
Fdv is piecewise polynomial. 

If X generates V, the Box Spline B{X) is a piecewise polynomial (relative 
to (X,A)) distribution supported on the zonotope Z{X). 



A REMARK ON THE CONVOLUTION WITH THE BOX SPLINE 3 

Let / be a smooth function on V. Then there is two distributions natu- 
rally associated to X, A, /: 

• the piecewise polynomial distribution B{X) *d f- 
on a test function test, 

{B{X) *d f, test) = / • • • / iest{X + ^ t^a^)dtl ■ ■ ■ cUn- 



• the smooth density B{X) *c /: 
on a test function test, 



c 



{B{X) *c f,test) = / f{v) / ••• / test{v + y2tiai)dti---dtNdv. 
Jv Jti=0 JtN=0 j^]^ 

The notations *d and *c means discrete, versus continuous. B{X) *d f is 
the convolution of B{X) with the discrete measure f{X)5\, while B{X)* 
f is the usual convolution of B{X) with the smooth density f{v)dv. The 
subscript *c is just for emphasis. The operation *rf is denoted *' in [5], [6] 
and is called semi-discrete convolution. 

Our aim is to write an explicit formula for the difference B{X) *d f — 
B{X) *,/. 

We also associate to a € X three operators: 

• the partial differential operator 

{daf){v) = jj{v + ea), 

• the difference operator 

{Vaf){v) = f{v)-f{v-a), 

• the integral operator 

{Iaf){v) = C f{v-ta)dt. 



Jo 

The operator is the convolution -B([a]) *c / with the Box Spline asso- 
ciated to the sequence with a single element a. 

These three operators respects the space of polynomial functions C[V] on 
V. The Taylor series formula implies that, on the space C[y], the operator 
la is the invertible operator given by 

-, „g oo -. 



In particular, if / € C[V] is a polynomial 

- e 

'da 



(2) BiX)*,f=lill^-^^)f]dv. 
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If /, J are subsequences of X, we define the operators dj = Hae/ 
Vj = rifeeJ They are defined on distributions. 

Recall that dYB{X) = VyB{X \ Y), if y is a subsequence of X. 

A subsequence y of X will be called long if the sequence X \ 1" do not 
generate the vector space V. A long subsequence Y, minimal along the long 
subsequences, is also called a cocircuit: then Y = X \ H where H is an 
admissible hyperplane. 

In our formula, when / is a polynomial, B(X) *d f — B(X) *c f is nat- 
urally expressed in function of the derivatives dyf with respect to long 
subsequences Y. 

2. PlECEWISE SMOOTH DISTRIBUTIONS 

Our aim is to write an explicit formula for the difference of the two distri- 
butions B{X) *a f and B{X) *c /■ As the first one is a piecewise polynomial 
distribution, the second a smooth density, we will need to introduce an 
intermediate space of distributions. We will use "piecewise smooth distri- 
butions" . Let us give a definition. 

We continue to assume that X generates V. 

Definition 2.1. A generalized function 6 on F will be called piecewise 

smooth (relative to X, A) if: 

• the generalized function b is locally L^, 

• on each tope r, there exists a smooth function 6(r) on the full space 
V sTich that the restriction of 6 to r coincides with the restriction of the 
smooth function 6(r) to r. 

In this definition, given a tope r, the function b restricted to r (as well as 
all its derivatives) extends continuously to the closure of r. However, these 
extensions do not always coincide on intersections of the closures of topes. 

If b is piecewise smooth, we then say that the distribution B := b(v)dv 
(given by integration against the locally function b) is piecewise smooth. 

It is clear that if we multiply a piecewise polynomial distribution B by 
a smooth function, we obtain a piecewise smooth distribution. Remark 
that the space of piecewise smooth distributions is stable by the operators 
Va, and by convolution with Box Splines B(Y) iY any subsequence of X). 
However, it is not stable under operators da- For example, dxB{X) = 
Vx-B(0) is a linear combination of 5 distributions. 

3. Multiple Bernoulli periodic polynomials 

Let U be the dual vector space to V and F C ?7 be the dual lattice to A. 
If y is a subsequence of X, we define 

UregiY) = {u G [/ ; (a, u) 0, for all a G y} 

and 

J- reg (y) = Fn[/reg(y). 
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Consider the periodic function on V given by the (oscillatory) sum 



(3) W{X){v) = Yl 



This is well defined as a generalized function on V. In the sense of gen- 
eralized functions, we have 



(4) dxWiX)iv) = Yl 



e 



7errcg(x) 

We will use this equation to construct "primitives" of parts of the Poisson 
formula. 

We will call the series W{X) a multiple Bernoulli series. Multiple Bernoulli 
series have been extensively studied by A. Szenes [7j. They are natural gen- 
eralizations of Bernoulli series: for A = Zcj and X^, := [w, w, . . . , w], where 
u} is repeated k times with k > 0, the series 

„2iiTnt 

is equal to —^B{k, t—[t]) where B{k, t) denotes the k^^ Bernoulli polynomial 
in variable t. In particular, for A; = 1, we have W{Xi){tuj) = ^ —t + [t] (see 
Figure [2]). 

We recall the following proposition f7] (see also [2], [1]). 

Proposition 3.1. IfX generates V , the generalized function W{X) is piece- 
wise polynomial (relative to {X,A)). 

Thus we will also call VF(X) a multiple periodic Bernoulli polynomial. 

The above proposition is proved by reduction to the one variable case. 
Indeed, the function yi — —/ — r can be decomposed in a sum of functions 

ILaeX \°-'^/ 

with respect to a basis a,-, of V extracted from X. This reduces 

the computation to the one dimensional case. A. Szenes [7] gave an efficient 
multidimensional explicit residue formula to compute W{X). 
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Example 3.2. Let V = Mei © Me2 with lattice A = Zei © Ze2. Let X = 
[ei, e2, ei + 62]. We write G F as u = v\ei + ^262. 
We compute the generahzed function 

^2iTr{nivi+n2V2) 

W{vi,V2) - Yl (2i7rni)(2i7rn2)(2i7r(ni + 77,2)) ' 

Then is a locally L^-function on V, periodic with respect to Zei +Ze2- 
To describe it, it is sufficient to write the formulae of W{vi,V2) for < 
vi < 1 and < ^2 < 1, which we compute (for example using the relation 
1 ^ 1 I 1 ) as • 

nin2 (ni+n2) m(ni+n2)^ 712(11+712)^ ' 



T^(V1,U2) = 



-|(1 + t"! - 2i;2)(t^l - 1 + t^2)(2ui - ^2), t"! < V2 

-g(ui - 2i;2)('yi — 1 + V2){2vi - 1 — ■U2), ^^1 > V2. 
Thus we see that is a piecewise polynomial function. 



Remark 3.3. If X do not generate V, W{X) is not locally L^: take X = 0, 
then, by Poisson formula, 1^(0) is the delta distribution of the lattice A. 

Definition 3.4. A subspace s of F generated by a subsequence of elements 
of X is called X-admissible. We denote by TZ the set of X-admissible sub- 
spaces of V. We denote by TZ' the set of proper X-admissible subspaces. 

The spaces s = 1/ and s = {0} are among the admissible subspaces of V. 
The set TZ' consists of all admissible subspaces of V, except s = 1/. 

Let s be an admissible subspace of V. Let us consider the list X\s, where 
we have removed from the list X all elements belonging to s. The projection 
of the list X \ s on l//s will be denoted by X/s. The image of the lattice A 
in y/s is a lattice in V/s. If X generates V, X/s generates V/s. Using the 
projection V V/s, we identify the piecewise polynomial function W{X/s) 
on V/s to a piecewise polynomial function on V constant along the affine 
spaces V + s. 

Define ;7rcg(X/s) = U^egiX \ s) n s-^. Thus Tj-^giX/s) :=Tn C/reg(^/s) is 
the set of elements 7 G P such that: 

(7, s) = for all s G s; (7, a) / for all a G X \ s. 

Identifying the dual space to V/s to the space S"*-, we see that the function 
W{X/s) is the function on V given by the series (convergent in the sense of 
generalized functions) 

Eg2j7r(i),7) 
.er..(x/s)n..x\s2-(a,7) 

This function is periodic with respect to the lattice A, piecewise polynomial 
on V (relative to X,A) and constant along v + s. 
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If s = V, the function W{X/s) is identically equal to 1, while if s = {0}, 
we obtain back our series W{X). 

4. A Formula 

Let us now state our formula. We assume, as before, that X generates V. 
For each s G 7^, we consider all possible decompositions of the list X \ s 
in disjoint lists / U J. If / is a smooth function, the function 

F{v) = W{X/s){v){diVjf){v) 

is a piecewise smooth function on V. If y is a subsequence of X, the 
convolution B{Y) * Fdv is well defined and the result is a piecewise smooth 
distribution on V that we denote by B{Y) *c {W{X/s)diV jf). 

Theorem 4.1. Let f he a smooth function on V . We have 

B{X)*df-B{X)*^f =Y, {-l)\'\B{{Xr\s)Ul)*,{W{X/s)diVjf). 

sen' icx\s 

In this formula J is the complement of the sequence I in X \s. 

This equality holds in the space of piecewise (relative to {X,A)) smooth 
distributions on V, relative to {X, A) . 

Remark 4.2. If / is a polynomial, the term B{X) *c/ is a polynomial density 
and all terms of the difference formula are locally polynomial distributions 
on V. 

Before proceeding, let us comment on the proof. As in [3] (see also [5]), 
we use the Poisson formula to compute B{X) *d f. Then we group the 
terms in the dual lattice F in strata according to the hyperplane arrange- 
ment Uagx{ct = 0}. We then use the Bernoulli series as primitives of the 
corresponding sums. This way, we introduce the needed derivatives of the 
function /. 

Proof. Let TZ be the set of admissible subspaces of V. We have the disjoint 
decomposition: 

(5) [/= □ C/reg(A/s). 

sen 

Let test be a test function on V. We compute 




We apply Poisson formula to the compactly supported smooth function 

q{w)=f{w) / B{X){v)test{w + v) 
Jv 

as our sum S is equal to J2xeA '?(-^)- obtain 

S = y [ e^'^^'"'^U{w)dw. 
-yerJy 
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The lattice F is a disjoint union of the sets ri.cg(^/s) = F n C/reg(-'^/s) 
associated to the admissible subspaces s. Remark that the set associated to 
s = y is {7 = 0}. The term in S corresponding to 7 = is Jy q{w)dw, that 
is {B{X) *e f,test). 

As in the generalized function sense 

7er,eg(x/s) 

we obtain 

S = J2 f WiX/s){w)i-l)\''\^\dx\sqiw)dw. 

The function q{w) is product of the two smooth functions f{w) and Jy B{X){v)test{w+ 
v). By Leibniz rule, 

5 := ^2 [ [ W{X/s){w)dif{w)B{X){v)djtest{w+v)dw. 

We first integrate in v and use the equation satisfied by the Box Spline 
{B{X),dbh) = -{B{X \ {b}),V.bh). 

Thus we obtain 

S = y] V (-1)1^1 / / W{X/s){w)dif{w)B{X\J){v){\7-jtest){w+v)dw. 
s iuj=x\s ■'^■'^ 

Let us integrate in w. We use the invariance of the integral by V;, : 
Jvi^hfi){w)f2{w)dw = Jy fi{w){S/^bf2)iw)dw. As 6 € J is in A, and 
W{X/s){w) is periodic, 

S = Y^ (-1)'^'/ [ B{X\J){v)W{X/s){w)Vjdif{w)test{w+v)dw. 

Writing TZ = {V} UTZ' , we obtain the formula of the theorem. 

□ 



On the space of polynomials, one has 

Vj5// = (H ^—^)dx\sf 
beJ 

if / U J = X \ s. 

Recall that the space D(X) of Dahmen-Micchelli polynomials is the space 
of polynomials on V such that dYf = for all long subsequences Y. In 
particular, if s is a proper subspace, the sequence X\s is a long subsequence. 
So if /, J are such that / UJ = X\ sand/e D{X), then Vjdjf = 0. 

As a corollary of our formula, if p S D{X), we see that B{X) *d p = 
B{X) *cP- Let us state more precisely this result of Dahmen-Micchelli [4] 
(see also [6], chapter 16). 
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Corollary 4.3. IfpeD{X), then 

P{v) := B{X) *dP = Y.p{X)B{X){v - A) 
A 

is a polynomial function on V , equal to (HaGX ^~%~~)P — ^i-^) *cP- 

In this formula, we have identified B{X), B(X) *dP, and B{X) *cP to 
piecewise polynomial functions. 

5. Vertices of the arrangement and semi-discrete 

convolutions. 

We now give a twisted version of Theorem 14.11 where we twist / by an 

exponential function _ 

The set of characters on A is the torus T := U /T . If 5 € T, we denote by 
the corresponding character on A. More precisely if g has representative 

G eU, then by definition g^ = e^'-^'^^'^K Define 

X{g) ■.= {aeX;g^ = 1}. 

U g £ T = U /r has representative G € f7, we denote hy g + T the set 
G + T. 

For a £ X, introduce the operator 

{Via,g)f){v) = f{v)-g-'^fiv-a). 
If y is a subsequence of X, define 

Vf,= J] V(a,5). 

aGV 

We introduce a subset TZ{g) of admissible subspaces, depending on g. 

Definition 5.1. The admissible space s is in TZ{g) if the space + T) n s-*- 
is non empty 

Remark that if G is not in F, then V is not in the set TZ{g). 

Remark 5.2. If s € 7^(5), then all elements of X n s are in X{g). Thus TZ{g) 
is contained in the set of admissible spaces for X[g). However the converse 
does not hold: take F = Mw X = [2a;], A = Zw, and G = \uj* . Then 
X[g) = X , so that V is an admissible subspace for X{g). However, V is not 
in n{g). 

If s G Tl{g), take G (5 + T) n s^. Then (5 + F) n s-^ is the translate by 
5s of the lattice F n s^. 
Define 

Freg(X/s,5) = {g + T)f^ C/reg(X/s)^. 
Thus Freg(Ar/s, g) consists of elements ^ G 5 + F such that 

{i, s) = for all s G s; a) 7^ for all a G X \ s. 



10 MICHELE VERGNE 

The following series 

is well defined as a generalized function on V. 

The function W{X/s, g){v) is not periodic with respect to A. We have 
instead the covariance formula 

(7) W{X/s, g){v - A) = g-^W{X/s, g){v). 
In the sense of generalized functions, we have 

(8) dx\sW{X/s,g){v)= Yl 

We recall the following proposition (see also [2],[1]). 

Proposition 5.3. The generalized function W{X/s,g) is a piecewise poly- 
nomial (relative to {X,A)) function on V. 

It is proven similarly by reduction to one variable. 

Example 5.4. Let V = Moj, A = and Xk := [uj,uj, . . . ,uj], where uj is 
repeated k times with A; > 0. Then F = Zw*. Then if z is not an integer 

WiXk,zu;*){toj) = Y 



We have, for example, (see [2]) 



^2iTT[t]z 



W{X,,zun{tu)- 

t-[t] , 1 



W{X2,zu*){tu) = e2-W^(- SL + 



Here [t] is the integral part of t. This function [t] is a constant on each 
interval ]£, i + 1[, and W{Xk, zco*) is a locally polynomial function of t. 

Theorem 5.5. Let G G U , and g its image in U/r. Letf{v) = e2*'^<^'<^>/i(z;), 
where h is a smooth function. Then 

B{X)*df= Y 5] (-i)i'ii?((xns)u/)*e(Ty(x/s,<7)5/V5^). 

senig) icx\s 

In this formula, J is the complement of I in X \ s. 

Remark 5.6. If G G F, then B{X) *d f = B{X) *d h, and the formula of the 
theorem above coincide with the formula of Theorem 14. II for h: the set TZ{g) 
coincide with the set TZ, and the term corresponding to V in the formula of 
Theorem [53] is B{X) *c h. 
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Proof. We proceed in the same way than the proof of Theorem 14.11 Let 
test be a test function on V. We compute S := Jy{B{X) *d f){v)test{v) by 
Poisson formula. If 



we obtain 



q{w) = h{w) / B{X){v)test{w + v), 
JV 



7Gr ' 



Thus 



CG{9+r) 

The set (7 + r is a disjoint union over s G Tl-ig) of the sets Trf.^{X/s,g) = 
{g + T)r\Urcg{X/s), so 

S= I W{X/s,g){w){-l)\''\^\dx\M{w)dw. 

Then, using Leibniz rule for da , and equation for the Box Spline, we obtain 
that S is equal to 

E E (-1)''' / / W{X\s,g){w)dif{w)B{X\J){v){V.jtest){w+v)dw. 

s£Tl{g) IUj=X\s •'^ •'^ 

Using the covariance formula ([7]) for \ s, (7), we see that 



W{X\s,g){w)h{w){V^bf2){w)dw= / W{X\s,g){w){V{h,g)h){w)f2{w)dw 
V JV 

and we obtain the formula of the theorem. □ 



Let us point out a corollary of this formula. 

Definition 5.7. We say that a point g G U /T is a toric vertex of the 
arrangement X, if X{g) generates V . We denote by V{X) the set of toric 
vertices of the arrangement X. 

If 5 is a vertex, there is a basis cr of ^ extracted from X such that 5" = 1, 
for all a £ a. We thus see that the set V{X) is finite. 
If X is unimodular, then V{X) is reduced to g = 0. 

Corollary 5.8. (Dahmen-Micchelli) 

Let g G V{X) he a toric vertex of the arrangement X and letp G D[X{g)) 
he a polynomial in the Dahmen-Micchelli space for X{g). Assume that g ^ 0- 
Let /(A) = g^p{X). Then B{X) *d f = 
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Proof. We apply the formula of Theorem 15 . 51 with h = p. As g ^ 0, all terms 
s € TZig) are proper subspaces of V. Let us show that all the terms in our 
formula are 0. Indeed let lUJ = X\s. Let /' = InX{g) and J' = JnX{g). 
Then I' U J' = X{g) \ s is a long subset of X{g). As Vf, = V//, we see that 
S/'Vj, is already equal to 0. 

□ 

I wish to thank Michel Duflo for comments on this text. 
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